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Recently it has been shown that cumulants significantly simplify the analysis of multipartite weak 
measurements. Here we consider the mathematical structure that underlies this, and find that it 
can be formulated in terms of what we call the moment algebra. Apart from resulting in simpler 
proofs, the flexibility of this structure allows generalizations of the original results to a number of 
weak measurement scenarios, including one where the weakly interacting pointers reach thermal 
equilibrium with the probed system. 



I. INTRODUCTION 



For many readers, the word "cumulant", if it means anything, probably evokes a slight feehng of discomfort: a 
, recollection, perhaps, of a baffling definition and a paper left half-read. Yet, as we hope to show here, cumulants 
should have pleasurable associations. They arise as part of an algebraic structure, the moment algebra, that can be 
defined very simply yet has striking properties. It has the familiar operations of complex analysis - a multiplication 
■ and inverse, functions like log and exp, a derivative operation, power series expansions, etc. - but all transposed into 
a very different setting, with functions defined on a lattice of finite subsets instead of the continuum of a complex 
space, and with curious-looking new definitions for the functions. In the moment algebra, the cumulant is just the 
C ' log function, though many textbooks do an impeccable job of concealing this fact. 

Cumulants have a long history, with roots in statistics. They were probably first considered by Thorwald N. Thiele 
[2] as "half-invariants" ; they then went through a protean sequence of name changes [1, 0, Q until the current name 
Q finally stuck. Since they are tools of statistics and probability theory , it is perhaps unsurprising that they have 
been applied in statistical mechanics [1, [^, [^, [l3| , notably for the calculation of virial coefficients and perturbation 
I expansions of the free energy, as well as in solid state physics, quantum chemistry, and quantum field theory (see e.g., 
^ i [ill, [l^. [13I. [1^ ) . Other studies have focused more directly on the cumulants as, in some sense, genuine multipartite 
correlation measures, both in a classical setting [l^, and a quantum setting p^ . 
, The second component of this paper is weak measurement fvA [isj . This is a way of obtaining information about 
CO ' a system while perturbing it only a little, by coupling the system weakly to a pointer. The imprecision of the 
, measurement outcome is compensated for by running many repeats of the protocol, each time with a freshly prepared 
I ' system. The measurement results obtained this way are often surprising ^T^, and give insight into the underlying 
00 \ physics, as in the analysis ^20| of Hardy's paradox [2l| . Other examples include phenomena in fiber optics (2^ and 
. photonic crystals psj . 

Here we consider multipartite weak measurements, by which we mean measurements involving more than one 
pointer [l^j HI, H^, H^l • The moments of pointer observables (i.e. the expectations of products of those observables) 
' turn out to depend in an extremely complicated way on weak values. Despite this, it has recently been shown that 
^ I the cumulants of pointer moments are very simply related to cumulants of weak values psj . This suggests that the 
■ - - ' role cumulants play in simplifying perturbation expansions in statistical mechanics may have an analogue in weak 
measurement. Note that, given the weak value cumulants, one may if one wishes obtain the weak values themselves 
by the inverse operation to the cumulant (the exponential in the moment algebra); this gives an operational procedure 
for computing weak values. 

With the help of the moment algebra, we uncover some of the mathematical structure that underlies the favourable 
interplay between cumulants and weak measurements. We also show that these results can be generalized by relaxing 
some of the assumptions behind weak measurement. For instance, we consider a situation where weak measurements 
are performed over an extended period during which the system undergoes continuous evolution. This is related, by 
what can be broadly described as an imaginary time transformation, to another scenario where the pointers and the 
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system correlate via thermalization; we call this "thermal weak measurement" . Thus in general we get a broader view 
of weak measurement, together with some tools for making its analysis more tractable. 

II. THE MOMENT ALGEBRA AND CUMULANTS 

Here we introduce what we call the moment algebra. This can be regarded as a natural setting for discussions 
about cumulants, and provides a handy formalism for the proofs in the rest of this paper. 

Let fin be the set of integers 1,2, ... ,n, for some n > 0. Let 7W„ denote the set of functions that assign a complex 
number /(a) G C to every subset a of r2„, including the empty set 0. We will refer to such an / as an Al-map. 

If / and g are two A^-maps in A4n, we can define their product fg simply as the new A4-map [fg]{a) = f(a)g{a), 
for any subset a of f2„. However, there is another product, the convolution product f * g, which has particularly 
interesting properties. If a = {oi, . . . , Ok}, let daf denote the formal derivative {d/d^ai ■ ■ ■ 9/9^0^)/, where the are 
notional variables that we never deal with explicitly. We now define 

if*g)ia)^daifg), (1) 

where the right hand side is interpreted as follows: In the standard expression for the derivative of the product fg, 
we make the replacements dbf f{b) and dbg g(b) for any subset b C a; i.e. we replace the derivative dbf by the 
value /(&) of the yVI-map f on b, and similarly with dbg. We also replace plain / by /(0). As an example, suppose 
a = {1,2}. Then we have 

diAfg) = (5i,2/).9 + {dif){d2g) + {d^m.g) + fid^^g), (2) 
and after replacing the derivatives we obtain the convolution product evaluated at a = {1,2}, 

(/ * g){l, 2) = /(I, 2)5(0) + /(1)5(2) + /(2)g(l) + /(0)g(l, 2). (3) 
In general, ([T]) gives the explicit rule 

(/*ff)(«)= E (4) 

aiUa2— a 

for any subset a ^ 0, where the sum runs over all ordered bipartitions (ai, a2) of a, including (0, a) and (a, 0) (treated 
as distinct). In the case a = we find (/ * 5)(0) = /(0)ff(0)- The vector space of A^-maps, A(„, together with the 
convolution product becomes an algebra, which we also denote by M.n and call the moment algebra. 

The terms "moment" and "convolution" are inspired by the following construction. Let (/>(a;i, . . .Xn) be a complex 
valued integrablc function of n complex variables. Define an A^-map by 



f4,{a)= J Xak---Xa^(f>{xi,...,Xn)dxi...dXn, (5) 

for each subset a = {ak, . .. ,ai}. Hence, the A^-map assigns to each subset a the moment of the function 
corresponding to a. If we now choose another function ip, then * gives moments of the usual convolution 
(j) * ipiui, ■ ■ ■ , Vn) = / 4>{xi,. . . , Xn)'>Piyi - xi, . . . ,yn - Xn)dxi . . . dxn- Wc Can write this succinctly as f^=¥ ^ f^^^, 
where the subscript '0 * ?/;' is the usual convolution product. In fact, any A4-map can be represented by an in 
this way, for some (non-unique) 0, but we do not make use of this, and proceed entirely within our abstract algebra 
framework. 

There is an identity, 1* in A(„, defined by 



r(0) = l, l*{a)=0 for a^$. (6) 



We have 



/*l* = l**/ = /. (7) 
Given a sufRciently differentiable mapping : C ^ C we can define a mapping F* : A^„ — > A^„ via 

(F*/)(a)-a.F(/), (8) 
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where we assume that the formal derivative operates according to the standard chain rule, giving daF{f) as a function 
of the formal derivatives dtf of/. Note that F*f{^l) — F(/(0)). Equation ([5]) defines an operation taking us from 
a function F : C ^ C to a function F* : 7W„ ^ Mn- Furthermore is a homomorphism under composition; viz., 

{FGy^F*G*, (9) 

for any two functions F, G : C C. 

We can readily extend the definition ^ of the "*" operation to operate on functions F with several complex 
variables, e.g., if F : — > C, and we have two A^-maps / and g, then 

F*if,g)ia)=d,Fif,g). (10) 

If we now choose G(/, g) = fg, we can use the extended definition to find G*{f,g) = f * g. We can furthermore apply 
dSl) to obtain 

F*if.g) = {Fifg)r. (11) 

Thus enables us to carry over maps of complex numbers to Mm preserving all their properties. For instance, if 
F{f) — f"^, we obtain an inverse in A4n, defined whenever /(0) by 

(/-!*)(«) = a,(ri), (12) 

where is the multiplicative inverse l//(a). Thus / * f^^* = 1*. Two other operations are 

(log*/)(a)=5,log(/), (13) 
(exp*/)(a)=5<,exp(/), (14) 

where "log" here and elsewhere means logg. From ^ we deduce that exp*(log*) = 1*, and from (fTT|) that 

log*(/*5)-log*/ + log*5, (15) 

and log* (/-!*) = - log* /-^ 

The prescription for the inverse, (|12p . can easily be turned into an explicit rule by formal differentiation. One finds 



/-i*(0) - l//(0), (16) 
ri*(l) = -/(l)//(0)^ (17) 
/-i*(l, 2) = -/(1, 2)//(0)2 + 2/(l)/(2)//(0)3, (18) 



and in general 



peA(a) •' 

the sum being taken over all partitions A(a) of a, with \p\ denoting the number of parts (subsets) of the partition p. 
This follows at once from the combinatorial version of Faa di Bruno's rule [2^ (and the moment algebra acquires an 
added grace by this use of the theorem of a beatified mathematician). It is remarkable that the simple definition in 
(HI) leads to such a rich structure. 
Similarly, Faa di Bruno's rule gives 

log*/(a)= ^ (b|-l)!(-l)l''|-i5£|^, a^0, (20) 

which is well-defined for A^-maps / such that /(0) 0. From (|13p and ([^0)1 we find 

log*/(0)=log/(0), (21) 
log*/(l) = /(l)//(0), (22) 
log* /(I, 2) = /(I, 2)//(0) - /(l)/(2)//(0)2. (23) 
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We can make a similar calculation for exp* and find 

exp* /(a) - e/(0) ^cepf{c)- (24) 

peX(a) 

We shall refer to log* / as the cumulant of the A^-map /, and exp* / as its anticumulant. Cumulants were originally 
introduced in the context of statistics and probability theory [?[ . To formulate these "classical" cumulants within this 
framework we let Xi, . . . ,X„ be a collection of random variables, and define an 7W-map / by /(a) = (Yij^aXj), for 
each subset a C {1, . . . , rt} = ri„, where (• • • ) is the expectation value of the product of the random variables. (We 
also take /(0) = (1) = 1.) Then 

log*/(a)=log*(n,eaX,) (25) 

is precisely the classical joint cumulant of the random variables {Xj}j^a, and exp* /(a) is their classical anticu- 
mulant. This justifies us in using the same terminology for log* / and exp* / in the more general situation where / 
is an arbitrary A^-map. 

We now use the moment algebra to rederive some properties of cumulants in this more general setting. Suppose 
that an A^-map / satisfies 

/(c) = /(cnA)/(cnB), VcCl]„, (26) 

where {A, B} is a bipartition of i.e., A n i? = and AiJB = 0„. We say that / factorizes with respect to {A, B} 
on ri„. Then log* /(c) = for any nonempty c C i^n such that c A ot c <^ B, and in particular log* f{^n) = 0. In 
other words the cumulant of factorizing A^-maps vanishes. To see this, we note that we can write / = /a * /s, where 
/a(c) = /(c) if c C A and /a(c) — otherwise, and /b(c) = f{c) ii c C B and /b(c) — otherwise. Then 

log* /(c) = log* fA{c) + log* /b(c), (27) 

and this is zero for any c that is not cither entirely in A or entirely in _B, since then there is some element a of c with 
a ^ A, and shows that log* /^ must vanish; but this is also true of log* Jb- It can be shown that this property 
characterizes cumulants [30, _31j . 

We mention here a property of cumulants that we will frequently make use of. Given a £ C, define the scalar a* by 
a*(0) = a, a*(c) = if c 7^ 0. Thus in the special case a = 1, a* is what we have previously called 1*, so our useage 
is consistent. Now / * a*(c) — {af){c). Thus the convolution product with a* is equivalent to scaling the value of / 
for all subsets of ri„ by the factor a. Note also that log* a*{c) = 0, for all c ^ 0. Thus, by ITS)) . 

log*(/ * a*)(c) = log* a*(c) + log* /(c) = log* /(c), (28) 

for c 7^ 0, so scaling / by a constant factor leaves the cumulant unchanged on all non-empty subsets of r2„. 

The comparison between A^„ and complex analysis is further strengthened by the existence of the power series 

log*(l*-H/)^/-^+ ^*^*^ ~..., (29) 
which converges whenever |/(0)| < 1. When applied to the empty set this becomes 

log*(l* + /)(0) = /(0)-^ + ^-..., (30) 

which is the familiar power series for the complex function log(l + /(0)), as it should be according to (HI]). When 
(I29p is applied to sets other than 0, checking its validity is a pleasant exercise. For instance, using the product rule 
(HI to evaluate the repeated convolutions gives 

log*(l* + /)(1, 2) = /(I, 2) - [/(I, 2)/(0) + /(l)/(2)] + [/(I, 2)/(0)2 + 2/(l)/(2)/(0)] + . . . 
_ /(1, 2) /(l)/(2) 
l + /(0) (l + /(0))2' 

which one sees is the correct expression if one compares it with (|23p and bears in mind the definition of 1* by ([6]). 
As a final ingredient, we define 

{d*f){a)=da{d,f),\iiia. (31) 

(The annoying restriction i ^ a can be removed by means of the multiset formalism in Section IVIIII ) This operation 
behaves like a partial derivative; for instance, we have 9* log* /(a) — (5*/) * f^^*{a). (In the case of the Al-map 
defined by (O, di corresponds to the operation (p — > d(j)/dxi.) It can also be interpreted as a sort of "raising" operator, 
taking /(a) to /(a U i). This lets us immediately derive log*(l* + f){a) for any a from log*(l* + /)(0), which, as we 
have just seen, is the complex function expansion (|30p . More generally, we can derive a moment algebra equivalent 
from any complex function power series. 
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III. WEAK MEASUREMENTS 

Weak measurement, developed by Aharonov and his colleagues, (TtI. [T8|. is a strategy for extracting information 
from a quantum system S, and has a number of distinct ideas behind it: 

(i) A measuring system V is weakly coupled to the original system S, and V is then uncoupled and measured. This 
allows some limited information about S to be gained with little disturbance to S. 

(ii) By repeating the entire procedure many times, with the system identically prepared on each occasion, the noisy 
information about S obtained by the weak coupling of V can be averaged to give a definite answer. 

(iii) The system S can be both preselected and postselected, the latter meaning that the procedure concludes with a 
measurement on 5, and only if a particular outcome is obtained are the data included in the averaging process. 

These ideas combine constructively and enable one to probe a system in new ways. In particular, one can express 
a weak measurement result in terms of a quantity called the weak value, akin to the standard expectation value. 
This depends upon both the pre- and post-conditioning states, and can take very curious-seeming values, that can 
nevertheless be shown to have a natural physical meaning [l8'|. 

In the standard weak measurement setup the pointer system V has a continuous degree of freedom, like the position 
of a single particle. Another common assumption is that this pointer particle initially is in a pure state (j) for which the 
expectation value of both the position and momentum observables are zero. Hence, the probability density for 
finding the particle at position a; is a gaussian centered at zero. Furthermore it is assumed that the coupling between 
the pointer and system is given by the impulsive Hamiltonian H — jp (g) AS{t), where p is the momentum operator 
on V and A is any Hermitian operator on S, and where d{t) is the delta distribution centered at time t = 0. Suppose 
the system is initially in state jV'i) and is postselected in state \ipf). The state of the pointer after the interaction and 
postselection is proportional to exp(—ijp (S) A)\(f))\'ipi). After the interaction, the pointer position q is measured, 
and by averaging over many repeats one obtains the expectation value to any desired accuracy. 

To model the weakness of the interaction we expand the resulting expectation value up to the first order in the 
interaction parameter 7, resulting in [l^ 

(g) = -fReA^, (32) 
where A^ is the weak value of the observable A defined by 

This basic type of weak measurement can readily be generalised. The pointer can be an arbitrary quantum system. 
We do not necessarily have to have a continuous degree of freedom. The Hilbert space could be finite-dimensional; 
e.g., the pointer could be the spin degree of freedom of a particle. The coupling can he H = jS{t)s (8) A, where s is 
now any Hermitian operator, and likewise one can measure any Hermitian operator r on V. We further allow the 
initial pointer state \(p) to be arbitrary. Then ([32)) becomes [1^ 



(r) ^ {r),p+jRe{(A^), (34) 

where 

C = -2z((rs)^-(r)^(s)^) (35) 

and (r)0 = {(j)\r\(j)). We can conclude that, by measuring the expectation value of the observable r on the pointer, we 
can obtain the weak value {A)^ on the system. 

A different direction of generalisation is to weakly measure several operators Ak either simultaneously 0, [2^ [2^ 
or sequentially [131 ■ In the latter case, one couples pointers at successive times tk via the Hamiltonians Hk = 
lk^{t — tk)sk "X) Ak, and assumes that the system evolves by unitaries Uk between these times. One can then calculate 
the moment (ri • ■ -Tn), i.e. the expectation of the product of pointer measurements. It turns out that this can be 
expressed in terms of sequential weak values given by 

, ^ {4'f\Un+lAnUn.--AiUi\,p,) 
{lpf\Un+lUn ■ ■ ■ Ul\lpi) 

The expression for (ri • • - rn) in terms of sequential weak values is horrendously complicated [2^. It undergoes a 
striking simplification, however, if one looks at cumulants. We note that /(a) = (J^jQafj) is an A^-map, since it is 
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well defined for every subset a, if we add the assumption that /(0) = 1. We will denote the cumulant log* f(a) by 
We can also define another A4-map using sequential weak values on subsets as 



A ( , {i>f\Un+lFnUn - ■ ■ F2U2FlUl\4^,) / A, if J G 

'"^"^^ (^^|C/„+iC/„...t/i|V..> ' "^^'"li ifj^a. ^-"'^ 

If we compare this with p6p we see that we only insert the operator Aj if j G a. As an explicit example, consider the 
case of four pointers and the sequential weak value ^^,(2,4), which would be 

We will refer to log* Ayj as the sequential weak value cumulant. 
The following theorem was proved in [28l |: 

Theorem III.l (Cumulant theorem). To the lowest joint order in the variables 7, 

log*(n,enr,) - (n,6o7j)i?e {^log* A^{n)} , (39) 

where ^ is given by 

^ = 2(-z)l"l (n,eo(r,s,)^^ - n,eo(r,)^^. . (40) 

As explained above, the parameters jj signify the interaction strength between the system and the pointers, and we 
obtain the weak measurement scenario by expanding the cumulant log* (Iljiz^rj) in the strength parameters 7, keeping 
only the lowest order coefficients. This is also the sense in which the equality in equation p9p is to be interpreted: 
as an equality up to the lowest orders. Note furthermore that the "joint order" of 7271 is 2, for 757372 it is 3, etc. 
Hence, in the above theorem the lowest (nonzero) joint order in the expansion is jfij = n. By this theorem we can 
conclude that the weak measurement setup can be used to measure the sequential weak value cumulant on the system, 
since in the weak limit the joint cumulant of the pointer observables is simply related to the sequential weak value 
cumulant on the system. This can thus be regarded as a natural generalization of the weak measurement scenario in 
the single pointer case, where the expectation value of the single pointer observable corresponds to the weak value on 
the system. 

In the case of simultaneous weak measurement (2^ . [25l . [26j , one couples n pointers at a time to via the Hamiltonian 
Hk = J2k ^(^ ~ ^o)7fc'Sfc ® A]., and the simultaneous weak values are given by 



1 ('0/l^7 r(n) ■ 



{A,,..., A,,U = ^ E ' 7; i"/^ ^ (41) 



where it runs over all permutations of 1, . . . , n. Again, one can define an A^-map and a cumulant denoted by log* A.u]s- 
There is an analogue to Theorem (jlll.ip for the simultaneous case [28]: 

Theorem III. 2 (Cumulant theorem for simultaneous measurement). // (Iljgo^'j) is the expected value of the product 
of n pointers in a simultaneous weak measurement, then, to the lowest joint order in the variables 7, 

log*(n,earj) = (n,eo7j)i?e {^log* A^^ift)} , (42) 

where ^ is as in Theorem UII.l]) . 



IV. A NEW CUMULANT THEOREM 



The original proof of Theorem IIII . 1 1 was by no means transparent. We will show how the moment algebra setting 
allows a better proof. However, we begin by considering a different way of gathering information from the pointers, 
where the corresponding theorem can be proved more simply (and in section rVIII we prove the original result, which 
requires an extra flourish of cumulant technology) . 

The assumption behind Theorem lIII.il is that the moment (Jlj^a^j) for a subset a of the total collection of pointers 
ri is obtained by coupling just that set of pointers to the system: in other words, to obtain (J\j^afj)^ one does an 
experiment in which just the pointers in a are coupled, so a separate experiment is needed for each subset. The 
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alternative approach that we now adopt is to suppose that a single experiment is carried out in which all the pointers 
are coupled and measured, but a subset of the measurement results is used for calculating each moment. 

Another way to put this is to say that in the present case we have the total time dependent Hamiltonian H{t) — 
1-p (8) Hg{t) + X]fe=i 7fe^(^ ~ tk)sk €5 Ak, where Hg(t) is the Hamiltonian that generates the unitary evolution between 
the coupling with the pointers. Hence, Hs{t) generates the sequence of unitary operators Uk- In Theorem IHI.H 
however, we had a separate experiment with a separate Hamiltonian Ha{t) = l-p ® Hs(t) + X^iea Ik^it ^ ii)si (8> Ai 
for each subset a. 

We shall see that the joint coupling of all pointers leads to a result that differs from Theorem lHI.il This comes 
about because of the perturbation of the system by those pointers whose readings are not being used to calculate 
(Hj-ga'^j)?;. When we measure a single pointer expectation {rk)ri this is given, up to first order, by the standard 
expression 7fci?e(^feA^), the coupling of other pointers producing perturbations that are only of second order in the 
7's. However, if we measure {rirj)^ in the presence of other pointers, the standard second order expression is no 
longer correct, other second order terms being introduced by other pointers. Despite these complications, we still 
obtain a succinct relationship between cumulants of pointers and weak values, as we shall see shortly. 

To state the theorem, let 77 be the state of all the pointers after coupling and the postselection on the system. 
(We define 77 more precisely in (|47l) below.) The moment of the pointer observables in the subset a is (Hjgarj),, = 
tr{riY\j^af'j), where the subscript 77 indicates that all the pointers are coupled. 

Theorem IV. 1 (Cumulant theorem with all pointers coupled). Suppose all n pointers are coupled sequentially. Let 
a he a subset of the finite collection of pointers il. To the lowest joint order in the variables 7, 

\og*{Ujear,)r, - iIljealj)Re {(log* A^a)} , (43) 

where ( is given by 

e - 2(-z)l-l l[i{r,s,)^^ - {r,)^A-'^)4',) = 2(-7)l'^l [] log*(r„ (44) 

jea jea 

This theorem allows us to consider the cumulants of subsets of the collection of pointers. This is a slight generaliza- 
tion of the theorems in Sec. IHIl and is more in line with the characterization in Sec. |lT]of cumulants as mappings from 
A^-maps to A^-maps. Just as the A^-map (Hjgarj)^ is defined for every subset of pointers a, so is the corresponding 
cumulant A^-map log*(Hjga''i}r)- 

Proof. First, we establish some notation. Let be the initial state of pointer k in state space Tik, and denote the 
total state space Tii ® . . . <8)7i„ of pointers by V, where we assume — n. Let be the initial state of the system 
S, and \^pf) the postselected final state. Define projectors on 5 ® by 

Pt = «) ([^^(^il • ■ • \(f>n){<l>n\) , 

Let rk be the pointer variable measured on pointer k and Hk = Ak be the coupling between pointer k and the 
system iS. For any state p on S ®V, define 

LkP=-i[Hk,p\, (45) 



Apart from the coupling to the pointer we also assume that the system S has a (possibly time dependent) Hamiltonian 
of its own, generating a unitary evolution between the couplings to the pointers. We let Uk be the unitary operator 
on S that gives the evolution between the coupling with pointer k — 1 and the coupling with pointer fc, and let Uq be 
the unitary that maps the initial state to the application of the first pointer. We then define 

UkP = (iv ® Uk)p{iv ® Uk)^ (46) 



to be the corresponding unitary channel. The state of the pointers after the interaction and postselection on the 
system can thus be written as 

^ tr5(P/^Y„+ieT"^"^„ ■ ■ ■ e-^^^^UiPj) 

tr(P/Z^„+ie7.i"Zi„---e7iiiZ^iP,) ■ ^ ^ 
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Let us define the A^-map (Iljgarj)^ = tr [{Ilj^arj)r]]. Suppose now that we replaced one of the observables rk in 
(Ilj(:a'rj)ri with the identity operator 1^. We would then obtain a new Ai-map Uk(a) = (Rj£a\{k}'''j)ri- This A^-map 
factorizes on a with respect to the partition ({A;}, a\{fc}), i.e., for all & C a we have Uk{h) = Uk{bC\{a\{k}))uk{hf^{k}). 
Thus we can conclude that log* Uk{a) =0. If we combine this with the multilinearity of the cumulant in the pointer 



observables we find that 



Define 



log*(njea(fj - (''i>0,))r, = log*(nje„rj)^. (48) 



|%|t/„+i.-.[/i|^.)P ■ ^^^^ 

The denominators of (|49p and (|T7)l do not depend on a, and are thus scalars in Ain- The scale-invariance of the 
cumulant, (I28p . thus yields 

log*(n,e,(r, - (r,)^J), = log* v. (50) 

We shall now expand the above cumulant to the lowest joint order in the parameters jj. Let us write for the 
derivative with respect to the variables "fj with labels in the set 6, i.e., 

(51) 

(These "proper" derivatives should not be confused with the formal derivatives da introduced in section |TT1) We 
next prove the following 

an * / w / if a \ & ^ 

log -(«)l7=o = |i,g*^(„) if (52) 

where 

\{lpf\Un+l- ■ ■Ul\tl)r)\- { 1 J^a, 

and 7 = means jj — 0, for 1 < j < n. Equation tells us that the first potentially nonzero expansion coefficient 
of the cumulant log* v{a) can itself be regarded as a cumulant, but of the new A4-map w. This method of regarding 
the expansion coefhcient of a cumulant as a cumulant in its own right is a technique that we will use repeatedly. 

To prove ([5^ . let us first suppose that a\b ^ 0. Hence, there must be some element j G a such that j ^ h. Recall the 
expression for the cumulant in terms of partitions, (j20p . and suppose p is a partition of a. Then 9^ncgpw(c)|-),=o — Oi 
since for that c G p that contains j the derivative d/d^j is not applied to v(c), and consequently this term contains 
rj — {rj)lj but not Lj and must therefore vanish. This gives the first part of ([5^ . For the case b = a, we again use 
the fact that, for any j ^ a, a term containing rj — {rj)lj but not Lj must vanish, to find that d2^ci£p'v{c)\y^o — 
Ilci£pdcv(c)\j^o, for any partition p of a. The statement in (|52p follows. 

It remains to evaluate the cumulant log* w{a) . Let us rephrase the definition of Lj in (j45p as 

L,=Lf + Lf'^ Lf'{p)^{~^s,^A,)p, Lf'\p) = p{ts,®Aj), (54) 
so the subscript 'left' ('right') indicates which side the operator is applied to. We can write 

^(") = Yl ^ci,c2(a), (55) 



Cl,C2 



where the sum is over all ordered bipartitions (ci ,02) of a, and where ,C2 (o) is defined as in (|53p but with Lj replaced 
by Lj°^' when j S ci, and by L"®''* when j G C2. By using the fact that and L"*''* act independently, and that Pi 
is a projector onto pure product states, one can show that Wci,c2 factorizes on a with respect to the partition {ci, C2}. 
Thus log* 

^ci,c2(^) — except when either ci — or C2 — 0; so log w[a) — log Wa,0(^) ^' log ''^0,a(^)- Direct 
calculation shows that log* j = (— i)l"l{nfega((?'fcSA;) ~ (7'fc)(sfc))}^^(a), and log* ^ gives the complex conjugate. 
The theorem follows. □ 
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V. SIMULTANEOUS WEAK MEASUREMENT WITH SYSTEM EVOLUTION 



The last section focussed on sequential weak measurement. In the case of simultaneous weak measurement, it 
is assumed [13, [H, HI] that the Hamiltonian has the form Hk = J2k ~ to)lkSk ® -^k, and the evolution of the 
system, which occurred between coupling of pointers in sequential weak measurement, can be ignored here since the 
coupling occurs impulsively and simultaneously for all pointers, and any evolution before or after the coupling can be 
incorporated into the initial and final states by writing Itpi) — Ui\ipi), (ipf \ = ('0/|C^2- 

But suppose the coupling occurs over a finite time. There is essentially no change in the analysis if we assume a 
Hamiltonian Hk = f{t) X]fe=i 9kSk ® ^fc, where f{t) defines some time-course for the coupling. However, there is now 
the possibility of having the system evolve while the coupling is occurring. Suppose this evolution is given by the 
Hamiltonian Hg. Then the total Hamiltonian is H — l-p ^ Hg + f{t) fffcSfc ® Ak- We assume here that f{t) is 

constant on a time interval of length r, and find it convenient to let our expansion parameters 7^ be the total strength 
of the interaction between the system and pointers, given by 7^ = J gkfit)dt = rg^. For times < i < r we can thus 
write the Hamiltonian as 



H=lv®Hs + yik—®Ak. (56) 



k=l 



As before we assume the system and the pointers start in a total product state, and that after the time interval r 
we postselect the system in \ipf)- The state of the pointers after the postselection is 



and we measure the observables ri, . . . ,r„, respectively, on each pointer, and calculate the corresponding cumulants 
log* {Iljf=arj) a- Now, for a = {oi, . . . let 

(A,, , . . . , J4r,+i , . . . , n] = (V-.le-^^IV'.) ^^^^ 

and define the A4-map 

^ „ k+l 

^(")^~feX! / (Aa,(,,, . . . ,^a,(i,)t«[Tfe+i, . . .ri](5(r - ^T,-)dTi • • -drfedrfe+i. (59) 

V can be regarded as an average over all possible sequential weak values of the operators Aa-^ , • ■ • , Aa,. , where we 
take all possible rearrangements of the order in which these operators are measured, as well as varying the time steps 
between the applications. 

Theorem V.l (Weak simultaneous measurement with system evolution). Let a he a subset of the finite collection of 
pointers i7. To the lowest joint order in the variables 7, 

log*(H,ear,), = (H,6,7,)i?e {^log* V{a)} , (60) 

where ^ is as in Theorem \IV. 1[ 

To prove this, we begin with the following 
Lemma V. 2. Let a — {ai, . . . ,ak}- Then 



^ JTfc + i,Tfc,...,Ti>0 

(61) 

where the sum is over all permutations tt of the set {1, . . . , A:}. 

Proof. We use the Dyson series (proof: differentiate both sides with respect to t) 

^riY+V) ^ ^rY ^ T ^ir-t.)Yy^t,Y ^ T f'' e'^r-t,)Y y ^it,-U)Y y ^t,Y ^^^^^^ (g2) 

Jo Jo Jo 
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Putting V = Y^^^iljXj, the only term that survives the combined operations of differentiation by and setting 
7 = is the /c-times repeated integral 



9X(^+^)|,=o = E f r ■■■ r e^^~'''^^X,.,,A-''-^''^X^ik^i)---^ (63) 
^ Jo Jo Jo 

To obtain (|6ip . make the change of variables ri = ti, r2 = t2 — ii, ■ • ■ , Tfc = tfc — ifc-i, Tk+i = t — tk- □ 

Proof of Theorem. Consider the A4-map log* {Hj^aTj) a- Following the proof of Theorem lIV.ll we can replace all 
the observables rj with rj — {rj)lj without changing the cumulant. Still following Theorem IIV.II we find that the 
lowest order term in the expansion has joint degree \a\ and corresponding expansion coefficient d2log* {Hj^aifj — 
{rj)ff,^))cr\j=o = log* d{a) with the new A^-map 

tr{Pf{n,eair, - {r,)i,))[dle-^-"]P.[dl^-"]}\^^^, (64) 



where the sum is over all ordered bipartitions (ci,C2) of a. Next, we apply Lemma rV.2l to both [d^-^e '"^^j |-y=o and 
j-g-y gjrffj 1^^^^ Y = -ii-p (g) Hs,Xk = -^Sk «) Ak, and Y = ii-p <S) Hs,Xk = ^Sk (X" Ak, respectively We find 
that d — D * D' , where 

D{a) - {-zp{U,^a{{r,s,)^^ - {r,)^^{s,)^,))V{a), (65) 

and where D' is the complex conjugate of D. The statement of Theorem IV. 1 1 follows from d ~ D * D' together with 
the fact that log* d = log* D + log* D' . □ 

The interpretation of this theorem is quite intuitive. If there were no evolution, i.e. = 0, then 2?(a) would be 
the simultaneous weak value, (2^. [25l. [26j and (j41[) . given by symmetrizing over all orders of applying operators, viz. 



^(") = ^ ' ■ ■ ■ ' ^-.(1))-' (66) 

where the factor 1/nl comes from integrating ri, . . . , r„, r„-)_i > with the constraint X)j=i ''j ~ When Hg is 
nonzero, we must in addition average over episodes of evolution under e^*^-^* between application of the Ak with the 
lengths of all episodes summing to r. 

The theorem implies that simultaneous weak measurement can be simulated by collections of sequential weak 
measurements, by sampling over permutations of the ordering of the applications of the pointers, as well as over the 
time steps between the applications of the pointers. 



VI. THERMAL WEAK MEASUREMENT 



In the previous section we stretched the concept of weak measurement a little by allowing the system to evolve 
while the pointers are coupled. Here we stretch it further by abandoning the notion of preselection and postselection 
(key ingredients of the original weak measurement philosophy [H, [s^l), and instead considering a system in thermal 
equilibrium. As we will see, this thermal weak measurement concept is, formally speaking, closely related to the 
simultaneous weak measurement with system evolution considered in the previous section. The correspondence 
between these two scenarios is analogous to that between path integrals and equilibrium systems under the imaginary 
time transformation t -i-^ it [l^ . 

For thermal equilibrium systems with Hamiltonian H, the Helmholz free energy [33 | can be written 

F = -ilogZ(/3) = -ilogtre-'5«, 

where /3 = 1/kT, with T being the temperature. When external parameters, e.g. fields, are changed infinitely slowly, 
the difference between the final and initial free energy is equal to the work performed on the system, under the 
assumption that the system is kept in contact with a heat bath at constant temperature T [3^ . The Taylor expansion 
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of the free energy with respect to the external fields thus characterizes the system's response to small changes. This 
picture can be extended to several fields gj coupling to the system via observables Aj , for instance with linear coupling 

Hc^Hs + Y^gjA,, (67) 
j 

leading to a free energy 

F = -^logtre-'^"''. (68) 

We follow (s^l and refer to the expansion coefficients of F with respect to 71, . . . ,7n, i.e. d2F\^^o for a subset a of 
the 7's, as generalized susceptibilities. We shall now construct a weak measurement scenario where the correlation of 
the pointers, as measured by the joint cumulant of the pointer observables, turns out to be directly proportional to 
these generalized susceptibilities. 

Instead of letting a collection of pointers weakly interact with the system for specific times, here we let the pointers 
and system equilibrate under the assumption of weak interactions. When this combined system has equilibrated we 
separate the pointers and, as before, measure a collection of observables on them. We therefore consider a total 
Hamiltonian 

n 

where 7^ = /Sgfc, and we assume that the system and the pointers reach the thermal equilibrium state under this 
Hamiltonian, yielding 

The expectation for pointer measurements is the A^-map {Ilj^a'''j)p = ^'''ip^jea'i'j)- We also define the Al-map 

g(a) - JI^^o (71) 

which gives, up to a normalising factor, the Taylor coefficients in the expansion of the partition function with respect 
to 7. 

Theorem VI. 1 (Weak measurement of a system in equilibrium). Let a be a subset of the finite collection of pointers 
fl. To the lowest joint order in 7 we find 

log*(n,ear,)p = (n,e,7,)eiog* £{a) - -l3{Iljeal,)id2F\^=o, (72) 

e = n,ea(tr(r,,s,) - tr(r,)tr(,s,)). (73) 

This assumes that the various traces tr(rj), tr(sj), and tr {rjSj) are well defined, which in the case of infinite 
dimensional Hilbert spaces requires them to be trace class The theorem tells us that the joint cumulant of the 
pointers is directly proportional to the generalized susceptibility of the system. 

As an application of these ideas, consider a collection of pointers in equilibrium with a heat bath. Correlations 
between the pointers will be generated by the heat bath, and these can be characterised by local observables. Our 
theorem says that, if the coupling of the pointers to the heat bath is weak, the cumulants of these local observables 
will be proportional to the generalized susceptibilities of the heat bath. 

As mentioned above, there is an analogy between thermal weak measurement and simultaneous weak measurement. 
To see this, we can use Lemma IV.2l to expand the A^-map £ in weak values as 

fc+i 

■dTk+i, (74) 



^^^^ = 4kYl 1 (^a„(,,,,--.,^a„(i))e[rfe+l,...,ri](S(/3-^rj)dTi--- 
P a "'rfe+i,...,ri>0 . 1 



where 



(Afe,...,^i)e[rfe+i,...,ri] = trie-l'^s] ' C^^) 

Comparing ([71)1 and (|75|) with (|59p and respectively, the t ^ it correspondence is clear; this allows us to carry 
over a large part of the proof of Theorem IV. II to the present theorem. 
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Proof of Theorem. Consider the A^-map {Iljfzarj)p = tr{pllj^arj). In the analogue of p9|) . instead of replacing rj by 
{rj — {rj)^,- Ij), we replace it by [rj — tr(rj)lj), the trace playing the role previously taken by the expectation. As before, 
this modification of the pointer observables does not change the cumulants, (Yij^arj) p — log* ijlj^aij'j —tr{rj)Vj) p. 
Defining 

njga(rj- -tr(rj)lj) 
^ - tre-PHs ^^^) 



we use scale- invariance, (l28l) . to show 

92(n,ea(r, - tr(r,)i)),|^.o = log* tr(e-^«Xp)|,.o (77) 

= log* 52tr(e-'^«Xp) 1^=0 (78) 
= log* tr(n,e,s,Xp) a2tr(e-''^<^)|^=o 
= .^log*f(a), 

where (1751) follows from ([77)1 by the same argument that derived ([5^ from Furthermore, 

eiog* £:(a) - Clog* a2tr(e-'^^<^)|^=o (79) 
= $9,^ log tr(e-''^'=) 1^.0 (80) 

where (|79p follows from scale-invariance, (pS)) . and ([50]) follows directly from the definition of log* using ([12]); see the 
Appendix for details. □ 

Note that in the proof of Theorem IV. 11 e"'"^^ operates on the left of Pi and its conjugate operates on the right, 
which leads to the real part, i?e{^log* X'(a)}, appearing in ((60l) . In the above theorem e"*^^ appears without its 
conjugate, so we get the whole of ^ log* £{a) in ([7^ . 

VII. NEW PROOF OF THE ORIGINAL THEOREM 

Finally, we shall show how our moment algebra methods can be used to prove the original theorem in 28]. Note 
that we here prove a slight generalization of Theorem lIII.il in the sense that we allow the cumulants to be taken over 
arbitrary subsets a of the total collection of pointers Vl. 

Proof of Theorem \III.1\ 

, . _ U{PfUn+iX.nUrr--XiUiP,) J r,e^^^^ if J £ a 

, . _ tr {PfUn+iYnUn---YiUiP^) ( ei=^= if 7 e a 

|(^/|[/„+i-..C/i|V.)P ' ^'"l i ifj^a ^^'^ 

z{a) = (n,6,r,) = (83) 

so the Al-map z{a) is the expectation for pointer measurements in the subset a. Note that the normalizations of 
the A^-maps x and y have been chosen so that x{%) = 1 and y{%) = 1. For y another convenient property is that 
y(a)l7=o = 1 for all a. 

The idea of the proof is as follows: if the ratio in (|83p were defined in terms of convolution operations, so we had 
z = X * y~^* instead of 2; = xy~^ , then this would imply log* z = log* x — log* y, leaving us with the much simpler 
task of calculating log* x and log* y. In fact, it turns out that, by expanding the Al-map z in powers of the 7's, we 
can achieve this switch from multiplicative to convolution operations: see (|86p. 

We first prove the equivalent of ([52]) for z: 

f)^}r^^*y(n\\ _/log*z(a)if 6 = a, . . 

d,\og z(a)|^=o- I if 6 c a, 5^ a ^^^^ 



13 



where z{a) = d2z{a)\^=o- To this end, for any b C a define the A4-map Zf, by zi,{c) = d'^^^z{c)\y^o for all c C a. By 

a,^log*z(a)|,=o = E (W-l)!(-l)l^l-'5,^n^Wl7=o 

= E (|p|-l)!(-l)l^l-in^&'nc<^)l7=o = log*z,(a), (85) 

pGTr(a) cGp 

using the fact that z{c) depends only on the 7's in c. Since Za ~ z, (18411 follows for the case b — a. lib ^ a, there is 
some j ^ a with j ^ b. When we put = 0, since ((85|) shows the term e^^^^ is not differentiated, the operator Lj 
does not appear in Zb{a). Therefore rj is not coupled to S and Zb factorises on a. We conclude that log* Zb{a) = 0. 

Equation tells us that log* 5(a) is the first non- vanishing term in the expansion of the cumulant log* {Tlj^aTj) 
in the 7's. In other words, the relevant expansion coefficient of the latter cumulant can be regarded as the cumulant 
of the new A^-map z. To calculate this new cumulant we use the usual law for differentiation of a product, together 
with the observation that if (61, 62) is bipartition of b then 



K^m,^o - (n,e^>.(r.)0,)[9,>(fei)]|^^„, 



7=0 



y(&2)^'^=°' 



which yields 



- y(b) 



7=0 biUb2=b 



62 



7=0 



where we have the new A^-maps 



= E (n,eb,(r,)^^)[9,X6i)] 

6iUb2=b 

hiUb2=b 
= 5 * j/"^*(6), 



x{b) = a^x(6)|^=o, 



92 
02 



y{h2) 



7=0 



V{b) = {^jeb{rj)4>,)d2y{b)\^=o. 
With the aid of we thus find that the cumulant of the 7W-map z can be decomposed as 

log* z{a) = log* x{a) - log* y(a). 
We next turn to the evaluation of log* x{a). Equations ([5T|) and (|57|) imply 

tr (PfUn+lXMn ■ ■ ■ XiUiP- 



x{a) 



\{^|:f\Un+l■■■UM^)? 



TjLj if J G a 
i if j ^ a 



(86) 
(87) 



(89) 



(90) 



This should be compared with ([55)1 in the proof of Theorem IIV.II Just as in that proof, where we defined 



Cl,C2) 

here we define ici,c2 where Lj is replaced by L^*^^* if j £ ci and by L"^*^' if j € C2. By the same argument, we find 
log* x{a) = log* Xafi{a) + log* :r0_a(a), and direct calculation gives 



Xa,%{a) = {-i)^°-\Tlj<za{rjSj)^^)Ayj{a), 
and i;0,a(a) gives the complex conjugate. Thus 

log*5;(a) =i?e{2(-i)l''l(nje,(r,Sj)^Jlog* A^,(a)}, 



(91) 



(92) 



and this gives (|39p with the first part of ^ (see PD| ). The evaluation of the cumulant of log* y{a) is analogous to the 
above, and results in the second half of ^. □ 
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VIII. MULTISETS AND MULTISET CUMULANTS 



Our definition of cumulants in Section was somewhat unconventional. A more standard definition for tlie classical 
cumulant is 

log*(n,ea^,) - 921og(e^-^^'='''=)l7=o, (93) 

which is readily seen (Appendix) to be equivalent to ours given by (j25p . Thus cumulants can be thought of as 
coefficients in a formal power series expansion in the 7's. To a statistician, the expression being expanded is the 
logarithm of the moment generating function. To a physicist, a comparison with the Helmholtz free energy, (j68p . is 
compellin g, a nd indeed there is a strong connection with thermodynamics [lol |. A combinatorialist can also stake a 
claim [36l.l37l|. 

The formal expansion for two variables begins 

log(,7iX,+,.x.^ ^ 71(^1) +72(^2) - - {X,f] - - -7172 [(^1^2) - (^i)(^2)] + . . . , (94) 

and the coefhcient of 7172 is familiar as the classical cumulant log*(A'iX2). However, this expansion also forces on 
one's attention terms in higher powers of the 7s, such as 7^((X^) — (Xi)^)/2. The moment algebra, as we introduced 
it in Sec. Ull does not include such terms. However, with a slight modification of the construction of the algebra, all 
of these higher order terms can be incorporated. The natural setting for it is not subsets of a set, but multisets [3^ . 
In a multiset, an element may occur any finite number of times. For instance, if the underlying set is = {1, 2, 3}, 
then {1,1}, and {1,1,1,2,3,3} are two examples of multisets in il^ (Note that, as with ordinary sets, the ordering 
does not matter, i.e., {1, 1, 1, 2, 3, 3} = {3, 1, 2, 3, 1, 1} .) 

To extend the moment algebra we define 7W-maps on multisets rather than subsets, so an A^-map / assigns a 
complex number f{a) to every multiset a of 17. The whole machinery of formal derivatives and their action on 
composite functions, as presented in Sec. [TTl goes through essentially unaltered. In particular, given an A^-map / 
defined on multisets, we can define higher order cumulants log* /(a) for any multiset a. 

As an example, consider the multiset that consists of {1}, {1,1}, {1,1,1}, etc.. Suppose / is a function on this 
multiset. Then we can apply ([8]) to a = {1, 1, 1}, for instance, so that from 

rl/l rll rl / rl\3 

5{M4}log/ = y-3^ + 2^, 
(with primes denoting differentiation by variable 1), we deduce 

log* f ({1 1 1}) - - smiilMiil) + 

As a special case, we define /({1, 1, . . . , 1}) ~ for some random variable X . Then we can rewrite the above 

V ' 

k 

expression as 

log* /({1, 1,1}) - {X^)-i{X^){X)+2{Xf^n:,{X). 



In this way we obtain the classical cumulant ^^{X), some others in the series being ki{X) = {X) and K2iX) 



{X'^) — {XY . These cumulants are widely used in statistics. They have nice properties; for instance, K,j{X) = for 



j > 3 is a necessary and sufficient condition for X to be Gaussian [39| . 

It is easy to see that any multiset cumulant can be obtained from our original cumulant with distinct variables simply 
by setting certain subsets of its variables equal. For instance, log* /({1, 1, 1}) can be derived from log* /({I, 2, 3}) 
by setting every '2' and '3' to a '1' (and leaving 'I's unchanged). Similarly log* /({1, 1, 2}) can be obtained from 
log* /({I, 2, 3}) by sending '3' to '1'. 

We can mimic this procedure in the case of multipartite weak measurements by treating subsets of pointers identi- 
cally; i.e., we can couple all the pointers in each subset via the same system observable and afterwards measure the 
same pointer observable on each of them. As a simple case, suppose we have just two pointers and we couple both to 
the system through the interaction Hamiltonian H = s ® A and finally carry out the same measurement r on them; 
note though we have to give labels to s and r to indicate which pointer they belong to; so we have si, S2 and ri, r2. 
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respectively, for pointers 1 and 2. To avoid the complication of ordering of the couplings of pointers, let us consider 
simultaneous weak measurement. Then Theorem IIIL2I gives 

(rirs) - (ri)(r2) - T^iie Uk2(A)„} (95) 

where we have suggestively written ^2(^)10 for the simultaneous weak value (^^)^ ~ (^-u;)^ given by lliTll . 

This gives us a procedure for gaining information about the multiset cumulants of weak values. With thermal 
weak measurement we can use this procedure to measure higher order susceptibilities. If we couple mi independent 
pointers to observable Ai, 1712 pointers to A2, etc, Theorem IVI.ll gives a direct relationship between the correlation 
of the pointers and the relevant susceptibility: 

IX. CONCLUSIONS 

In physics we often study the effects of weak coupling between systems. By focussing on the effect of one system 
upon the other, weak measurement gives a way of understanding the nature of such interactions. If we weakly couple 
two systems, S and V, sa^ and then carry out a strong measurement on V, the effects of the coupling can be expressed 
in terms of weak values [iTlllsj. When P consists of a very simple system, e.g., a "pointer" or particle in a given initial 
state, the measurement results depend on weak values in a very simple way ((33)) . As P becomes more complicated, the 
dependency becomes rapidly more complicated, and in fact already assumes a highly baroque form when V consists 
of two pointers applied at different times (see the Appendix of 28]). However, this complication vanishes if one takes 
cumulants: the cumulant of the measured variables and the cumulant of the weak values are once more simply related. 
The aim of this paper has been to give proofs of this fact that illuminate why this phenomenon occurs. 

The proofs of our various theorems repeatedly use two properties of cumulants: first, that they are logarithms in 
the moment algebra Mm and turn a convolution product into a sum; second, that they vanish on maps that factorise, 
i.e., that can be written as a product of two maps that are defined on disjoint sets of variables. All the maps that 
we construct are elements of which is therefore the natural setting for the proofs. This algebra is in itself an 
interesting object. Although some of its features have been thoroughly described in the literature, we are not aware of 
any explicit formulation of the algebra as an object in its own right. We feel that it deserves this recognition, because 
of the simplicity of the definitions of its operations, and the surprising richness of the structure that this gives rise to. 

Cumulants have long played a role in statistical mechanics [5l.[8l.[9l.[l0|. where they are used to simplify perturbation 
expansions. This suggests alternative weak measurement scenarios. In this spirit we consider a collection of pointers 
that reach thermal equilibrium with a probed system, and call this a "thermal weak measurement" . By its very 
nature, this excludes pre- and postselection, which are standard components of weak measurement. We lose thereby 
some of the strengths of weak measurement: many of the more intriguing phenomena in the standard setting arise 
from postselection. However, we retain the advantages of minimal perturbation of a system, and the possibility of 
applying several probes simultaneously or sequentially opens up some new territory for exploration. 

X. APPENDIX: DEFINING CUMULANTS VIA GENERATING FUNCTIONS 

Here we show that the standard definition of the classical cumulant via the generating function is equivalent to 
our definition (fT3| of log* / for the A^-map /(a) = (Hj^aXj). In other words, we wish to show that, for any multiset 
a, 

log*/(a)=521og/i(7)k=o, (97) 

where 

/i(7) = (e^"'"^')- (98) 

We first note that 



d2\ogh = A{h,djh,d]h,dl2,...), 



(99) 
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where A is a function of all the relevant partial derivatives of h, obtained via the chain rule when we apply to the 
logarithm. The same function A appears when we express the cumulant log* f{a) in terms of formal derivatives 

log* /(a) =9, log /(a) 

=A(/,9i/,92/,5i,2/,...) (100) 
=A(/(0),/(l),/(2),/(l,2),...). 

To obtain ^ from ^ and (fTOO]) we only need observe that d2h\^=o = (HjizcXj) = /(c). 

This equivalence of definitions can be extended to other functions h if we take /(a) — djhl^^o- For example, let 
h{j) = tr{e~^^^) = tr{e~^^^~^^ ^i^j)-)^ Then the above arguments prove (|80p . which occurs in the proof of Theorem 
IVI.ll Note that the arguments extend without difficulty to non-commutative observables. 
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